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ABSTRACT 

Rapid rotation in field red giant stars is a relatively rare but well-studied phenomenon; here we 
investigate the potential role of planet accretion in spinning up these stars. Using Zahn's theory of 
tidal friction and stellar evolution models, we compute the decay of a planet's orbit into its evolving 
host star and the resulting transfer of angular momentum into the stellar convective envelope. This 
experiment assesses the frequency of planet ingestion and rapid rotation on the red giant branch 
(RGB) for a sample of 99 known exoplanet host stars. We find that the known exoplanets are 
indeed capable of creating rapid rotators; however, the expected fraction due to planet ingestion is 
only ~10% of the total seen in surveys of present-day red giants. Of the planets ingested, we find 
that those with smaller initial semimajor axes are more likely to create rapid rotators because these 
planets are accreted when the stellar moment of inertia is smallest. We also find that many planets 
may be ingested prior to the RGB phase, contrary to the expectation that accretion would generally 
occur when the stellar radii expand significantly as giants. Finally, our models suggest that the rapid 
rotation signal from ingested planets is most likely to be seen on the lower RGB, which is also where 
alternative mechanisms for spin- up, e.g., angular momentum dredged up from the stellar core, do 
not operate. Thus, rapid rotators on the lower RGB are the best candidates to search for definitive 
evidence of systems that have experienced planet accretion. 
Subject headings: planetary systems — stars: evolution — stars: rotation 



1. INTRODUCTION 

The rotation rates of solitary field giant stars as a func- 
tion of temperature show a sharp transition from fast 
rotati on to slow rotation redward of the GOIII spectral 
type (|Gravlll9"89[ l. The majority of stars redward of the 
transition are characterized by a mean rotat ional veloc- 
ity of ~ 2 km s~^ ([de Medeiros et al.lll996l ). However, 
a small percentage of these red giant stars are found to 
have rapid rotation, which is generally defined as v sin i 
> 10 km s^^. Such high rotation rates are unexpected 
with normal stellar evolution because angular m omen- 
tum i s conserved as the stars expand. Moreover, iGravl 
(|1989( ) argues that the sharp transition from high to low 
rotation at GOIII indicates an ancillary braking mecha- 
nism caused by a magnetic dynamo effect created as the 
stellar convection zone deepens, further slowing the ro- 
tation; this mechanism renders rapid rotation even less 
likely in red giants. 

With the lack of obvious progenitors under normal evo- 
lutionary processes, these unusual, rapidly rotating red 
giants must be spun up by some mechanism, either in- 
ternal or external. One propose d internal mechanisrn 
is angular momentum dredge-up (jSimon fc Drakdll98l 
iFekel fc BalachandranI [T993f ). which relies on the star 
having a rapidly spinning core that rotates independent 
of the envelope. The core acts as an angular momen- 
tum reservoir that is tapped as the convective envelope 
extends deep into the star during first dredge-up on the 
red giant branch (RGB). 

An alternative mechanism is a gain in angular mo- 
mentum from an external source, which, in the ab- 
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sence of a cluster environment or stellar companion, 
would be the orbital angul ar momentum deliver ed by 
a brown dwarf or planet. iPeterson et al . (1983) were 
the first to consider planets as a source of angular mo- 
mentum in evolved s tars. Later studies (jSokeil Il998t 
ISoker fc Harpajl2000l ) invoked planet-induced mass loss 
as the possible "second parameter" determining hori- 
zontal branch ( HB) morphologies . In studies of red 
giant chemistry, lAlexandeil ()1967l ) first invoked planet 
engulfment by red giants to explain lithium enhance- 
ments in some of these stars. Since these early stud- 
ies, the idea of planet engulfment has been invoked re- 
peatedly explain both lithium enhance ment and rapid 

rotation observed in RGB stars (Siess fc Livid 
.. .. -—^f^ 



Reddv et al.ll2002t iDrake et al.l l2002t 1Carnev et all 



1999 



2003 



Denissenkov fc Her wigl l2004l and references therein) 



Siess fc Livi o (1999) modeled the accretion of sub-stellar 



companions by red giant stars and calculated observa- 
tional signatures of the accretion in those stars, such as 
infrared excess and ''Li enrichment. From the actual 
occurrence of these predicted observational signatures 
found in the red giant population, they estimated that 
4-8% of Sun-like stars must host planets if most planet- 
hosting (PH) stars have sig nificant interact i ons w ith their 
planets while on the RGB. lLivio fc Sokeil (|2002l ). on the 
other hand, used the extrasolar planets known at the 
time to estimate that at least 3.5% of red giants would 
have their evolution "significantly affected" by planets, 
with rotation-induced mass-loss being the most signifi- 
cant alteration. A comprehensive discussion of the effects 
planets may have on a star's evolution from the RGB 
through the HB and asymptotic giant branch (AG B) to 
the planetary nebula phase is given in [Soke? ('2OO I. 

In this study, we expand on the idea of Livio fc Sokeil 
(|2002f ) of using the known characteristics of planet- 



2 



Carlberg, Majewski, & Arras 



hosting stars to predict future interactions between the 
RGB stars and planets. We are particularly interested 
in trying to account for the frequency of RGB stars with 
rapid rotation by planets depositing angular momentum 
into their parent stars. Although wc focus on planets, 
we do note that brown dwarfs and low-mass stellar com- 
panions are also capable of spinning up giants. However, 
brown dwarfs turn out to be unlikely angular momen- 
tum sources because few stars have brown dwarfs or- 
biting near enough to their stars to tidally i nteract sig- 
nifica ntly with their host stars f|Grether fc Lineweaveii 
|2006( ). Low -mass stars are twice as co mmon as mas- 
sive planets (jGrether fc Lineweavei<200"6h : however, their 
higher mass affords them a higher survival rate in com- 
mon envelope evolution. Therefore, they may not explain 
rapid rotation in samples of field giant stars that were 
selected to be companionless. We return to the issue of 
low-mass stars in 

Using known properties of PH stars as a basis for ex- 
ploring giant rapid rotation is advantageous for a num- 
ber of reasons. First, it is not immediately obvious how 
some parameters affect the probability of a star becom- 
ing a future rapid rotator. As an example, consider the 
role of the semimajor axis of the planet's orbit. The or- 
bital radius must be sufficiently small that the planet is 
eventually engulfed by the star, but within this criterion 
one would naively suspect that planets with the largest 
possible orbital radius (and thus largest angular momen- 
tum for a given mass) would maximize the probability 
of spinning up the star. On the other hand, as we show 
here, planets orbiting at larger radii will be engulfed later 
in the evolution of the star, when the stellar moment of 
inertia is larger and subsequent spin-up is smaller. Gon- 
sequently, whether or not a PH star can become a rapid 
rotator depends on both the angular momentum avail- 
able in the planetary orbit as well as when in the course of 
the host star's evolution the planet is actually accreted. 
The relative values of these properties are specific to indi- 
vidual planetary systems; therefore, using the properties 
of known PH stars incorporates the natural distributions 
of both the stellar and planetary properties, all of which 
affect the census of rapid rotators. By combining the dis- 
tributions of PH stars and planets with stellar evolution 
models that allow us to track the individual evolution 
of the PH stars into the RGB phase, we can begin to 
ask questions about the relationship between PH stars 
and the current population of red giant rapid rotators. 
How many PH stars will become rapid rotators? What 
is the lifetime of a red giant rapid rotator that is created 
through planet accretion? Is the number and lifetime of 
rapid rotators from expected planet ingestion enough to 
account for the number of observed rapid rotators in the 
current red giant population? 

These are questions we address here. We begin by 
describing the properties of current red giant rapid rota- 
tors and PH stars as well as our adopted stellar evolution 
models in S}2l In fJSj we present our model for the orbital 
tidal evolution and ingestion of exoplanets; the results 
of our modeling are presented in ^ Finally, in fJ5] we 
discuss the possible effects of observational biases intro- 
duced by different exoplanet detection methods as well 
as the unmodeled effects of mass loss and magnetic dy- 
namos. 
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Fig. 1. — Stellar m ass and meta llicity for a sample of low-mass 
single red giants from IMassarotti~ et al. (2008). The different sized 
circles indicate bins of vs'mi < 4 km s~^ (small), 4 < vsini < 
8 km (medium), and i;sini > 8 km s~^ (large). The vertical 
line indicates the maximum stellar mass probed by our PH sample. 

2. STELLAR PROPERTIES 

2.1. Frequency of Red Giant Rapid Rotators 

The peculiarity of rapid rotation among the red gi- 
ant population has been well studied and documented 
by many authors. While the exact value of the rota- 
tional velocity separating "rapid rotators" and "normal 
rotators" is somewhat arbitrary, we follow the literature 
standard of defining a rapid rotator as any red giant with 
V sin i > 10 km s~^. In the catalog of ro tational veloci- 
ties from Gleb ocki fc Stawikowskil ()2000l ) . rapid rotators 
account for about 5 % of the red giant popula tion. In a 
sample of 900 stars, Ide Medeiros et all ()19960 find that 
less than 5% of the red giants have vsini > 5 km 
and faster rotators are even rarer. A recent study by 
J. Carlberg et al. (in preparation) of almost 1300 distant 
K giants found that less than 6% of their sample were 
rapid rotators (depending on the number of undetected 
binaries that masquerade as rapid rotators), while the 
st udy of almost 750 r iearby single giants from Hipparcos, 
bv'Massa rotti et al] ()2008f ). yields only 1% of giant stars 
with rapid rotation. The general consensus, as evidenced 
by the four studies mentioned, is that the frequency of 
rapid rotation in the red giant population is only a few 
percent. 

In this work, wc compare predictions of the frequency 
of rapid rotation from planet ingestion (as well as trends 
of rapid rotation frequency with stellar properties) with 
the sa mple of observed red giants in the lMassarotti et al.l 
((200l survey because the Hipparcos distances allow 
accurate estimates of absolute stellar luminosities and 
masses. Figure [T] shows the masses and metallicities of 
669 stars from this survey grouped into three velocity 
bins. Included are all stars from the survey with masses 
between 0.7 and 10 Mq, temperatures less than ~5550 K, 
and < 0.001 (as suggested by the authors of that pa- 
per to avoid possible undetected binaries). In this sam- 
ple, ten stars (1.5%) show rapid rotation, 122 (18.2%) 
are in an intermediate v sin i bin (arbitrarily defined as 
4 < vsini < 10 km s~^), and the remaining stars are 
slow rotators. For reference, the vertical line shows that 
maximum mass probed by the PH stars described in the 
next section. In the stellar mass range probed by our 
exoplanet host star sample, 5 stars (1.3%) show rapid 
rotation, 40 (10.7%) are in the intermediate vsini bin. 
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and the remaining stars are slow rotators. Note in Figure 
[1] that there is a paucity of both rapid and intermediate 
rotators for the most metal poor stars for the low-mass 
range to the left of the line. 

2.2. Exoplanet Data 

The properties of exoplanet host stars were down- 
loaded from The Extrasolar Planets Encyclopaedic?, 
which is a compilation of both the known orbital param- 
eters of exoplanets and the physical properties of their 
host stars. We selected all those exoplanet host stars 
for which masses and orbital radii of their planets were 
available, for which stellar masses were available, which 
were known to be luminosity class V stars and thus on 
the main sequence (MS), and which would evolve to the 
RGB in less than 13.5 Gyr (a total of 99 systems with 
115 planets). A second sample adds all of the stars with 
no known luminosity class but that meet the rest of our 
selection criteria (38 additional planets in 31 additional 
systems), and for these, we simply assume that they are 
likely to be MS stars. In the following paragraphs we will 
give statistics for the first sample followed by those of the 
second sample in parentheses. The PH stars have [Fe/H] 
ranging from -0.32 to 0.50 dex (-0.68 to 0.50) with an 
average of 0.09 dex (0.10), while the masses range from 
0.91 to 1.43 Mq (0.87 to 3.9) with an average of 1.1 Mq 

Eccentricities of the planetary orbits were used in our 
calculations when available, and we assumed circular or- 
bits in the absence of such a measurement. Because the 
majority of exoplanets have been discovered by the ra- 
dial velocity method, the observed planetary masses suf- 
fer from an inclination ambiguity. However, observable 
stellar rotational velocities suffer from the same ambi- 
guity. Therefore, in using the minimum planetary mass 
(MpSini) we can compute the observable rotational ve- 
locity, wsini, for the PH stars that ingest planets to 
directly compare to v sin i measurements of actual red 
giant stars. This substitution is valid under the assump- 
tion that the planets' orbital axes are aligned with their 
stars' rotational axes. For comparison, this assumption 
is valid in our own solar system because all of the planets' 
orbital axes are tipped by only few degrees with respect 
to the ecliptic, and the echptic itself is tipped by only 7° 
to the Sun's rotational axis. 

2.3. Evolution Models 

The evolution of a star is mostly determined by its 
mass and composition alone ^ . All of the PH stars in our 
sample have mass measurements, and all except 14 (18) 
have measured metallicities. We assume solar metallicity 
for those without such a measurement. With these two 
properties we assign each star a stellar evol ution track. 
We u se the stellar evolution models of Girar di et al.l 
(|2000f ). which provides a grid of models with metallicities 
ranging from Z = 0.0004 to Z = 0.03 ([Fe/H] « -1.74 to 
0.18 dex) and a stellar mass range of 0.15-7.0 M©. The 

^ http://www.cxoplanet.eu/catalog.php, database for "All Can- 
didates detected," accessed on 2008 July 17. 

^ Rotation also plays an important role; however, grids of evo- 
lution models explicitly accounting for rotational effects are not 
currently available for the range of stellar masses of the stars in 
our study. Some unmodeled effects of rotation are addressed in 
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models give values of the stellar age, effective tempera- 
ture and luminosity (hence radiu s), evolution stage, and 
qdisc as a function of time. (See lGirardi et al.]|2000l for 
a complete description of the models and all available 
global parameters.) 

The parameter qdisc is the mass fraction contained 
within the radius where the chemical composition first 
differs from that at the surface. We use it as an indica- 
tor of the amount of mass in the homogenized convective 
envelope, given by M^nv = M*(l — qdisc)- The mass 
within the convective envelope is important because this 
is the regio n of the star where tidal dissipation occurs 
()Zahnlll977i ): it affects the rate at which a planet spirals 
in and is the region where the angular momentum is de- 
posited. The qdisc parameter works well as an indicator 
of mass in the convective envelope from the MS turn-off 
up to first dredge-up on the RGB because the convection 
zone is responsible for keeping the star homogenized from 
the surface to the discontinuity, qdisc- However, once 
first dredge-up ends the base of the convection zone re- 
cedes back towards the surface, while the star remains 
chemically homogenized to the deepest point previously 
reached by the convection zone. At this evolutionary 
stage, qdisc no longer traces the convection zone. It re- 
mains static until the outward-moving hydrogen burning 
shell reaches this point, after which qdisc follows the H- 
burning shell. Thus, we find that qdisc becomes an imper- 
fect tracer of the convective envelope after first dredge- 
up and at that point gives an overestimate of the mass 
in the convective envelope. This leads to an underesti- 
mate of the stellar v sin i after first dredge-up for stars 
that accrete planets. We also find that qdisc is an im- 
perfect tracer of the convective envelope on the main 
sequence, where the mass in the convective envelope is 
generally on the order of a few percent of the stellar mass 
while qdisc would give a value between 20-90%. For each 
stellar mass, we use a constant value for the fraction of 
mass in the convective envelope when the star is on the 
main sequence, which we take from the flat regions of the 
plots of stella r enve lope mass versus age in Figure 1 of 
iMurrav et aP (|2001h . 

3. EVOLUTION OF STELLAR ROTATIONAL VELOCITY 

3.1. Angular Momentum Budget in Planetary Orbits 

The observed rotational velocity of the PH stars will 
evolve as the moment of inertia of the stellar envelope 
evolves or if the total angular momentum in the con- 
vective envelope changes. The moment of inertia in a 
stellar convective envelope, assuming an n = 3/2 poly- 
trope model, is given by lenv = (l/8)Me„i,i?^ , where 
Menv is the envelope mass and i?* is the stellar radius. 
The moment of inertia will steadily increase as the star 
ascends the RGB, primarily from the increase in i?*. For 
simplicity, we assume that the total angular momentum 
in the envelope can only change if a planet is accreted; 
we ignore mass lost during the stellar evolution, which 
is also ignored in the stellar evolution models. When a 
planet is accreted, the change in the stellar rotational ve- 
locity, Awrot, from the deposition of the planet's orbital 
angular momentum can be calculated as 



8[Afpv/GAf,ap(l-e^)] 
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The terms and factors in Equation ([T]) include the plan- 
etary mass (Mp), semimajor axis (flp), and eccentric- 
ity (e), the stellar mass (M^,), convective envelope mass 
(Menv), and radius (i?*), and the gravitational constant 
(G). The expression in square brackets in Equation ([T]) 
is the initial orbital angular momentum of the planet. 

Planets sufficiently close to the star will experience a 
strong hydrodynamical drag force in the red giant's en- 
velope, leading to orbital decay. Tidal friction with the 
red giant's envelope causes orbital decay for larger orbits 
where hydrodynamical effects are negligible. All plan- 
ets with s emima j or ax es less than a radius, ri„, will be 
accreted. ISoked (|1998| ) reports that a good approxima- 
tion for determining which planets will spiral in during 
a star's ascent up the RGB is Tin 5Rup, where Rup is 
the stellar radius at the tip of the red giant branch. This 
conclusio n was based on tidal equations originally pre- 
sented in ISokeij ()1996[ ) and is valid for stars with masses 
between 0.8 and 2.2 Mq. For a solar mass star that 
reaches 200 Rq at the RGB tip, all planets within ~5 
AU constitute an angular momentum reservoir that will 
be accreted during the RGB lifetime based on Soker's 
approximation. Applying this approximation to the PH 
stars, we can calculate the angular momentum and Avrot 
gained from accreted planets with Op < 5 AU. In Figure 
[21 we plot this Avrot, calculated when the PH star is 
both at the base of the RGB and at the tip of the RGB. 
Points for these two evolutionary stages for the same PH 
star are connected by a vertical line. On average, an ad- 
ditional 8 km in rotational velocity, indicated with 
the horizontal line, would put planet-ingesting giant stars 
(which generally rotate at ~2 km s~^ before ingestion) 
into the rapid rotator regime. This analysis illustrates 
two important facts. First, the angular momentum in 
planetary orbits is sufficient to create RGB rapid rota- 
tors in most cases. Second, the time of planet accretion 
with respect to the stellar evolutionary phase is critical 
because the same planet that can create a rapid rotator 
if it is accreted at or near the RGB base will not cre- 
ate a rapid rotator if it is accreted when the star is at 
the RGB tip. Additionally, a star that does become a 
rapid rotator will "spin down" during its ascent up the 
RGB. Consequently, the point at which the planet is en- 
gulfed and the time that a star spends as a rapid rotator 
are very relevant to understanding how planet accretion 
contributes to the rapid rotator population. 

3.2. Orbital Decay by Tidal Friction 

To find when the PH stars will actually accrete their 
planets, we followed the tidal evolution of the individ- 
ual planetary orbital radii through the evolution of the 
PH star from the zero-age main sequence (ZAMS) to 
the RGB tip. Tidal interactions can cause a planet to 
spiral into the star if the rotation frequency of the star 
is less than the orbital frequency of the planet. (For 
the Sun, this condition would hold for planets orbiting 
within ^ 0.2 AU, whereas for a 15 Rq red giant with 
usini = 2 km s^^, this is true for planets out to 1 AU.) 
In such a scenario, tides raised on the star by the planet 
will lag behind the line connecting the centers of the two 
bodies, and the star will experience a tidal torque that 
will act in the direction that would increase the rota- 
tion frequency of the star. This spin-up drains angular 
momentum from the planet's orbit, and the planet con- 
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Fig. 2. — Additional rotational velocity that PH stars would 
gain on the red giant branch if the orbital angular momentum of 
all their respective planets with < 5 AU is deposited into the 
stellar convective envelopes. This is calculated both when the star 
is at the base of the RGB (triangles) and when it is at the tip of 
the RGB (circles). Solid vertical lines connect the results for the 
two evolution stages for the same star. Stars gaining more than 
8 km of added rotational velocity, indicated by the horizontal 
line, would be considered rapid rotators. In cases where connect- 
ing lines overlap (e.g., the three stars at 1.4 Mq), upper triangles 
connect to upper circles and lower triangles to lower circles. 

sequently moves closer to the star. The amount of lag 
and the strength of the torque depend on the rate that 
energy is dissipated in the tidally distorted body, and we 
adopt the m odel for turbulent dissipation presented in 
IZahnI (|1977f ). The time-scale for the inward spiral of the 
planet depends sensitively on the separation of the star 
and plai iet. We integrate th e tidal orbital dec ay eq ua- 
tions bv lVerbunt fc PhinnevI (|1995f) and|Zalin| l| 19771 ) to 
find the expression for the maximum orbital separation, 
r.im of accreted planets, which can be computed for every 
time step in the stellar evolution models. This separation 
is given by 

(2) 

The term I{t) is the integral of time-dependent stel- 
lar properties that affect the tidal evolution, calculated 
in years from the ZAMS {t' = 0) to the evolution- 
ary time step at t\ it is described by Equation (5) of 
IVerbunt fc PhinnevI ()1995f ) and reproduced here: 

(3) 

We compute the value of I(t) for every time step of the 
stellar evolution models through numerical integration. 

The exoplanets are accreted by the host star in our 
models when Op = rm- Because the tidal decay 
timescales in IZahnI ()1977() depend on {ap/R^)^, the 
planet spends very little time in the region between 
its original separation and the stellar surface; thus, the 
angular momentum is deposited almost instantaneously 
compared to the evolutionary timescale of the star. The 
planets will deposit 40-60% of their angular momentum 
simply by moving from their initial separation to R^ ; we 
assume the remaining angular momentum is added as 
the planet quickly evaporates in the envelope (see §5.11 
for a discussion of this assumption). In this experiment, 
the evolution is followed up to the first ascent on the red 
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Fig. 3. — Accretion radius, Tin, calculated using Equation (|2j 
for a Jupiter mass planet in a circular orbit around an evolving 
host star as a function of the stellar age. Solar metallicity stellar 
evolution tracks for stellar masses ranging between 0.8 A/q and 
1.6 Mq were used to follow the evolution from ZAMS to the RGB 
tip. Solid lines connect the individual time steps for each mass 
model, and evolution progresses from bottom left to upper right. 
The two gray points on each track denote special evolution stages. 
The yellow point is the exhaustion of central hydrogen and the red 
point is the base of the RGB. 

giant branch. All of the PH stars are of low enough mass 
that they go through the helium flash prior to landing on 
the HB. The helium flash is poorly understood and is not 
treated in evolution models so we limit our analysis to the 
better understood parts of stellar evolution. In Figure [3l 
we illustrate how rin increases during the stellar evolu- 
tion for solar metallicity stars with a Jupiter-mass planet 
in a circular orbit. We show the evolution from ZAMS 
to the tip of the RGB for stars with masses between 0.8 
and 1.6 M©. For stars of the same mass but different 
metallicity, the behavior of the Tin looks essentially the 
same. The only differences are the evolutionary time- 
scales (the lower metallicity stars evolve more quickly) 
and the maximum value of ri„, which is a factor of 2 
smaller for the most metal-poor model compared to the 
most metal rich. In FigureSl we show the RGB evolution 
of four of the stellar mass models in Figure [3] on a linear 
age scale to clarify trends for this phase of evolution that 
occurs more rapidly. 

4. RESULTS 

For the 99 PH stars of known luminosity class, we flnd 
that 89 stars accrete at least one planet, 10 stars accrete 
no planets, and 36 become rapid rotators at some point 
on the RGB by gaining Aw sini > 8 km s~^. The proper- 
ties of these 36 star systems are listed in Table [U which 
gives the name, mass and metallicity of the star, mass 
and semimajor axis of the accreted planet, the maximum 
rotational velocity the star attains on the RGB, the evo- 
lution stage of the star when the planet is accreted (ei- 
ther MS, sub-giant (SG), or RGB), and the fraction of 
the RGB lifetime that the star is a rapid rotator {Irr). 
For the two stars that have two known companion plan- 
ets, HD 168443 b and HD 202206 b are accreted in our 
models while HD 168443 c and HD 202206 c are not. 

In Figure El we compare and [Fe/H] of the planet- 
hosting stars that accrete planets, that do not accrete 
planets, and that are red giant rapid rotator progeni- 
tors (RRPs) to see whether any trends exist. In general, 
there are no strong correlations between M, or [Fe/H] 
and whether a PH star absorbs a planet. To better ex- 
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Fig. 4. — Expansion of the red giant branch evolution for four 
of the models presented in the previous figure. The 0.8, 1.0, 1.2, 
and 1.4 Mq tracks are shown. The large point indicates the base 
of the RGB. Note the different time-scales of each panel. 

plore the differences in these samples, we plot cumulative 
distribution functions (CDFs) in and [Fe/H] in Fig- 
ures |6] and [7l respectively. Plotted separately are CDFs 
of the entire PH stellar sample, PH stars that do not 
accrete planets, PH stars that accrete one or more plan- 
ets but are not RRPs, and the RRPs. The distribu- 
tions of stars that accrete planets diverge from the dis- 
tribution of all the stars in general in the mass range of 
1.05-1.25 Mq. In this mass range, a larger fraction of 
the stars that accrete planets do not become rapid rota- 
tors. The two-sided Kolmogorov-Smirnov (KS) statistic 
reveals that the probability that the RRPs and the stars 
that ingest planets but are not RRPs are drawn from 
the same parent population is only 4%. The stars that 
do not accrete planets, which is the smallest subsample 
comprises a larger fraction of low-mass stars than the 
parent sample. Comparing all stars that ingest planets 
to all those that do not, the KS probability of being from 
the same parent sample is 16%. 

The [Fe/H] CDFs in Figure □ show that the low [Fe/H] 
stars tend not to accrete planets as readily. For this 
plot, we removed all stars for which we had assumed 
solar [Fe/H]. Interestingly, this plot also shows that of 
the stars that do ingest planets, those that are RRPs 
are somewhat more metal poor than those that are not 
RRPs. Repeating the KS-test for the [Fe/H] distribu- 
tion, we find that there is a 54% chance that the RRPs 
and the stars that ingest planets but are not RRPs are 
from the same parent population. Stars that do not in- 
gest planets and stars that do ingest planets have a 19% 
chance of being drawn from the same population. The 
stars that do not ingest planets form the most metal poor 
subsample. A tentative conclusion is that both the A/» 
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TABLE 1 

EXOPLANET HOST STARS THAT ARE RAPID ROTATOR PROGENITORS 



Star Name 


M^, 


[Fc/H] 






Max. Atisinj 




Irr 




(Me) 


(dcx) 




(AU) 


(km s-i) 


(Stage) 


(%) 


70 Vir 


1.10 


-0.03 


7.44 


0.48 


15.4 


RGB 


2 


CoRoT-Exo-l'' 


0.95 


-0.30 


1.03 


0.03 


15.8 


MS 


14 


CoRoT-Exo-2'^ 


0.97 


O.OO'^ 


3.31 


0.03 


56.8 


MS 


70 


CoRoT-Exo-3'' 


1.27 


-0.05 


21.60 


0.05 


286.3 


MS 


98 


HAT-P-2'' 


1.30 


0.11 


8.62 


0.07 


115.1 


SG 


91 


HAT-P-S*^ 


1.16 


0.24 


1.06 


0.04 


14.3 


SG 


15 


HAT-P-6'' 


1.29 


-0.13 


1.06 


0.05 


14.2 


SG 


17 


HAT-P-9'' 


1.28 


0.12 


0.78 


0.05 


10.2 


SG 


16 


HD 136118 


1.24 


-0.07 


11.90 


2.30 


9.3 


RGB 


0.06 


HD 141937 


1.00 


0.16 


9.70 


1.52 


11.0 


RGB 


0.2 


HD 168443 


1.06 


0.03 


8.02 


0.30 


27.3 


RGB 


11 


HD 178911 B 


1.07 


0.28 


6.29 


0.32 


14.5 


RGB 


3 


HD 179949 


1.28 


0.22 


0.95 


0.05 


12.5 


SG 


16 


HD 202206 


1.13 


0.37 


17.40 


0.83 


28.0 


RGB 


2 


HD 209458*' 


1.01 


0.04 


0.69 


0.05 


9.9 


RGB 


12 


T T T~\ O O r* 


1.25 


-0.12 


9.10 


1.10 


9.3 


RGB 


0.1 


xlJJ Doyoo 


1 1 Q 

l.io 


n 0/1 
U.z4 


1 on 

i.yu 


U.U / 


1 Q 
io.O 


"D T3 


oo 
lA 




1.24 


0.29 


1.87 


0.04 


26.1 


SG 


34 


HD 86081 


1.21 


0.26 


1.50 


0.04 


20.7 


SG 


26 


OGLE-TR-IO'^ 


1.18 


0.12 


0.63 


0.04 


8.6 


SG 


15 


OGLE-TR-56'' 


1.17 


0.00"= 


1.29 


0.02 


12.0 


MS 


16 


OGLE-TR-132'' 


1.26 


0.37 


1.14 


0.03 


11.5 


MS 


16 


OGLE-TR-182'' 


1.14 


0.37 


1.01 


0.05 


17.5 


SG 


16 


SWEEPS-04'' 


1.24 


0.00^= 


3.80 


0.06 


63.1 


SG 


75 


Tau Boo 


1.30 


0.28 


3.90 


0.05 


51.5 


SG 


71 


TrES-2b 


0.98 


-0.15 


1.20 


0.04 


23.9 


MS 


36 


TrES-3'' 


0.92 


-0.19 


1.92 


0.02 


30.7 


SG 


41 


WASP-l'' 


1.24 


0.00<= 


0.89 


0.04 


12.1 


SG 


16 


WASP-3'' 


1.24 


O.OO'^ 


1.76 


0.03 


20.7 


MS 


26 


WASP-5'' 


0.97 


O.OO'^ 


1.58 


0.03 


27.1 


MS 


36 


WASP-14'' 


1.32 


O.OO'^ 


7.73 


0.04 


93.0 


MS 


88 


XO-l'^ 


1.00 


O.OO'^ 


0.90 


0.05 


20.3 


RGB 


27 


XO-2b 


0.98 


0.45 


0.57 


0.04 


11.0 


SG 


17 


XO-3'' 


1.21 


-0.18 


11.79 


0.05 


167.7 


MS 


93 


XO-4'' 


1.32 


-0.04 


1.72 


0.06 


25.5 


SG 


36 


XO-5'> 


1.00 


0.18 


1.15 


0.05 


25.0 


SG 


36 



In all cases, the accreted planet was the first discovered for that star, e.g., HD 168442 
b was accreted but HD 168443 c was not.'' Transiting system'^ No measured stellar 
[Fc/H] 




0.9 1.0 1.1 1.2 1.3 1.4 1.5 0.9 1.0 1.1 1.2 1.3 1.4 1.5 

M, (Mo) M, (Mo) 



Fig. 5. — Distribution of planet hosting stars that will not ac- FiG. 6. — Discrete CDF of the stellar masses of: all stars in the 

Crete any planets (squares), that will accrete at least one planet sample (solid), PH stars that will not accrete any planets (dotted), 

(triangles), and that will be spun- up enough by the accretion of PH stars that will accrete planets but not become rapid rotators 

one or more planets to become rapid rotators for part of their red (dashed), and PH stars that will be rapid rotators (dot-dashed), 
giant lifetime (circles). 

ets, they are more likely to be RRPs. Stars slightly more 
and [Fe/H] are correlated with the probability of a PH massive than the Sun that accrete planets are correlated 
star accreting a planet and becoming a rapid rotator, al- with a lower probability of becoming a rapid rotator, al- 
beit weakly. The metal poor PH stars are less likely to though it should be noted that this may be a reflection 
accrete their planets, but of those that do accrete plan- of the [Fe/H] effect; the stars in this mass range are on 
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' ■'^ I ' ' ' All siarT 

No accretion 

. Accretion (not RRP)-' 

Q_g Accretion (RRP) 





(no rapid rotation) 

Accreted 

(rapid ["otation) 



1.00 10.00 



100.00 



Fig. 7. — Discrete CDF of the metallicities of: all stars in the 
sample (solid), PH stars that will not accrete any planets (dotted), 
PH stars that will accrete planets but not become rapid rotators 
(dashed), and PH stars that will be rapid rotators (dot-dashed). 



Fig. 10. — Discrete CDF of the semimajor axes of: all the ex- 
oplanets (solid), exoplanets that will survive (dotted), exoplanets 
that will be accreted but not create rapid rotation (dashed), and 
exoplanets that will be accreted orbiting RRPs (dot-dashed). 



100.00 



10-00 r ®* ® U e ° -] 

I " I t i ' ABO T: - 

5. 1.00 . ® ill. , . ' nB 

' 'a ' ^ 

0-10: Planet not accreted □ E 

Planet accreted A 
' Max. RGB Ax)„, > 8.0 km s"' O " 

0.01 I , , , 

0.01 0.10 1.00 10.00 

Fig. 8. — Distribution of planet mass and semimajor axis of exo- 
planets that will survive (squares), exoplanets that will be accreted 
(triangles), and exoplanets orbiting a star that is a RRP (circles). 
IBecause some planets are in multiple systems, a planet that will 
survive might still orbit an RRP. 




0.01 0.10 1.00 10.00 100.00 

Fig. 9. — Discrete CDF of the planetary masses of: all the 
exoplanets (solid), exoplanets that will survive (dotted), exoplanets 
that will be accreted but not create rapid rotation (dashed), and 
exoplanets that will be accreted orbiting RRPs (dot-dashed). 

average more metal rich than the sample outside that 
range — 0.13 dex compared to 0.08 dex (excluding stars 
with [Fe/H] set to solar in both cases). 

The role of the planetary parameters in planet accre- 
tion and the creation of rapid rotators is illustrated in 
Figures [SVfTUl Figure [5] shows the distribution of Alp and 
Qp for the planets that survive the entire RGB evolution, 



planets that are accreted, and planets with host stars 
that are RRPs. Note that there are two planets that 
survive the RGB phase orbiting an RRP that is spun up 
by the accretion of a closer-orbiting planet. Most planets 
with Up < 2 AU are accreted. CDFs of all the planets, 
planets that survive, planets that are accreted but do not 
cause rapid rotation, and planets that do create rapid ro- 
tation are plotted for Mp (Figure [51) and Op (Figure [TU)). 
Figure [5] reveals that no planet with Mp < 0.5 Mj^p 
has enough angular momentum to create a rapid rota- 
tor. The CDFs in ap of these different samples in Figure 
[TU]show large disparity. As expected, planets that are in- 
gested generally have smaller semimajor axes than those 
that are not ingested. However, this figure also shows 
that the ingested planets that orbit RRPs have a distri- 
bution more heavily weighted to smaller semimajor axes 
compared to the sample of ingested planets that do not 
orbit RRPs, despite the fact that more distant planets 
have more angular momentum for a given planet mass. 
The KS test gives a 3% probability of these two samples 
being drawn from the same parent population. This re- 
sult is explained by the fact that the tidal decay timescale 
depends on (ap/i?*)*. The planet spirals into the star es- 
sentially when i?* reaches roughly some critical fraction 
of Op. Consequently, while the available orbital angular 
momentum is increasing with \/a,p, the total A^rot added 
to the star is decreasing with ^Ja^ (i.e., the behavior of 
Equation ([1]) for i?* cx ap). The spin- up effect from the 
larger angular momentum available in planets with larger 
orbital radii is more than offset by the increased stellar 
moment of inertia when these outer planets are finally 
accreted, and we see this effect in Figure [TUl 

In order to compare these results to the Galactic red 
giant population in general, we must have a sense of not 
only how many PH stars become rapid rotators but also 
how long they remain rapid rotators. Figure [11] shows 
the rapid rotation lifetime, Irr, as a fraction of the total 
RGB lifetime, for the 36 RRPs as a function of the model 
A/*; the approximate total RGB lifetime is given on the 
top axis. The average value of Irr is 31% ± 4%, and this 
is indicated by the horizontal line. Although there is a 
slight positive correlation between the stellar mass and 
Irr., there is also significant scatter in Irr, particularly at 
the highest and lowest masses (M* > 1.2Mq and < 
I.QMq). Thus, we assume that the average fractional 
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Fig. 11. — Fraction of the RGB lifetime that the planet-hosting 
rapid rotator progenitor will spend as a rapid rotator, lnn, as a 
function of the stellar mass. Approximate RGB lifetimes are given 
on the top axis. The dashed line indicates the mean value of Ijiji- 



TABLE 2 

Summary of results for subsets of the PH 

STARS. 





NpH (dp < 3 AU) 


Irrp 


Irr 


fexp 


Case A 


96 


37.5% 


31% 


0.58% 


Case B 


130 


33.8% 


31% 


0.53% 


Case C 


96 


33.3% 


33% 


0.55% 



lifetime is applicable to a generic population of red giants 
over a range of masses. 

The expected fraction of rapid rotators, fexp, pre- 
dicted from our planet ingestion models can be written as 
fexp = fpHfRRplRR, where fpH is the fraction of stars 
that host planets, fpRp is the fraction of planet-hosting 
stars that will become rapid rotators from planet inges- 
tion, and Ijiji is described above. While the total fraction 
of stars that harbor planets of any mass or at any orbital 
distance is unknown, it is possible to calculate the frac- 
tion of stars with planets within the observational con- 
straints imposed by exoplanet survey methods. In par- 
ticular, Grether & Lineweaver (2006) find that 5%±2% 
of all Sun-like stars have at least one planetary compan- 
ion orbiting within 3 AU. Using fpH — 0.05, we calculate 
fRRP only considering those stars in our sample with at 
least one planet within 3 AU. In our sample, the number 
of PH stars meeting this criterion, which we call NpH, 
is 96, and of these stars, 36 accrete a planet with suf- 
ficient orbital angular momentum early enough in their 
evolution to cause rapid rotation. Therefore, the fraction 
of RRPs, fpRp, in this sample is 37.5%. Assuming the 
stars spend an average of 31% of the RGB lifetime as 
rapid rotators, then the expected fraction of rapid rota- 
tors, fexp, among single red giants from the accretion of 
planetary companions is 0.58%±0.26%. The actual fre- 
quency of rapid rotation in various surveys of red giants 
ranges from 2-5%. As such, planet accretion can account 
for only 10-30% of these stars if the present known PH 
sample is a reasonable proxy for field RGB progenitors. 
These results {NpH, fpRP, Irr, and fexp) are tabulated 
in Tabled] under "Case A." 

If we include the stars with unknown luminosity class 
and assume that they are likely on the main sequence, 
we add an additional 34 stars with at least one planet 
orbiting within 3 AU, of which 31 ingest planets, and 8 



of these are RRPs. Following the same calculations de- 
scribed in the previous paragraph, we find 0.53%±0.23% 
of low-mass red giants should be rapidly rotating due 
to planet ingestion; these results are summarized under 
"Case B" in Table [21 

5. DISCUSSION 

5.1. Evaporating Giant Planets in the Stellar Envelope 

In our models, the issue of whether a planet is fully 
evaporated when it is accreted can affect the amount 
of angular momentum deposited by about a factor of 
2. As noted earlier, 50% of the planet's angular mo- 
mentum is added to the star simply from tidal decay of 
its orbit from the initial separation to the stellar sur- 
face. We assume the remaining angular momentum is 
contributed to the star through complete subsumption, 
and this assumption warrants further discussion. A com- 
plete physical treatment of the planet evaporation would 
result in a significant increase not only in the complex- 
ity of the physics in our models (particularly since we 
are modeling almost 100 individual systems) but also in 
the number of assumptions regarding unknown charac- 
teristics of these exoplanets, such as their physical radii, 
compositions, and density profiles. Instead, we look to 
previous mod els of the ev aporation of planets in main se- 
guence stars (iSandauist et al. 1998 . _2002) and giant stars 
(|Livio fc Sok eil ll984l:ISokedll998l) for some insight. 

For MS stars. ISandauist et al. ()1998[ ) modeled the ac- 
cretion of Jupiter- and Saturn-like planets into stars of 
1 Mq and 1.22 Mq and found that even though the two 
stellar models differed by only 0.22 Mq, the details of 
the stellar structure were different enough to change the 
amount of evaporation from total (1 M© case) to only 
one-third of the original planet mass (1.22 Mq case). 
Additionally, they found a sensitive dependence of the 
amount of evaporation on the structure of the plane- 
tary interior; planets with steeper density gradients had 
higher survival rates than those with sha ll ower density 
gradients. For giant stars, iLivio fc Sokeil (|1984f ) found 
that planets that spiraled into their stars could suffer one 
of three fates — complete evaporation, partial evapora- 
tion, or mass accretion. Their models covered a range of 
planet masses, initial separations, and physical assump- 
tions for stellar models at 0.88 and 1.2 Mq. They found 
that a critical planet mass, merit, exists above which a 
planet will accrete mass and below which a planet will 
evaporate. Depending on the model and physical as- 
sumptions, this critical mass was found to lie between 
~7.5 and 10 Mj„p. For the models that gave the smaller 
value, planets up to only ^--^5 Mj^p are evaporated com- 
pletely, while those with masses between 5 and 7.5 Mjup 
are partially evaporated. However, a later paper by 
iSiess fc Liviol (11999 1 gives merit 20 Mj„p. For stars 
near the RGB t ip ( which have the most rarified atmo- 
spheres) , [Soke3 ()1998l) found a relationship between the 
fate of planets and the planet mass. Planets of order a 
Jupiter mass or less evaporate at distances from the stel- 
lar center that depends on the planet mass. A Jupiter 
mass planet evaporates at approximately 10 Rq, and less 
massive planets evaporate at larger radii. Planets more 
massive than Jupiter may avoid evaporation until they 
are deep enough in the star to overflow their Roche lobes; 
therefore, the most massive destroyed planets may con- 
tribute angular momentum to the core rather than the 
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envelope. 

In light of all the uncertainties in these models, we 
conclude that our assumption of complete evaporation is 
not unreasonable, particularly for the lower mass planets. 
If we take a more conservative stance and assume that 
planets more massive than 5 Mjup do not evaporate, then 
these planets only add angular momentum to their stars 
that is equal to the difference between the orbital angular 
momentum when the planet is at and when the planet 
is at i?,. Applying this assumption to the PH stars we 
modeled in Case A, we find that four of the stars will no 
longer become rapid rotators; therefore, Jrrp drops to 
33%. Additionally, the average lifetime increases slightly 
to 33%, and fexp becomes 0.55% ± 0.25%. These results 
are listed under "Case C" in Table [2] 

5.2. Effects of Mass Loss on the RGB 

Ignoring mass loss will not affect the number of possi- 
ble giant rapid rotators if the stars accrete planets on the 
RGB itself; it only affects the time that star spends as a 
rapid rotator. However, a surprising result of our simula- 
tion is the number of planets for which ingestion occurs 
before the RGB. Figure [T2l shows on the Hertzsprung- 
Russell (HR) diagram the evolution of the PH stars that 
will become rapid rotators. The triangles mark the evo- 
lution stage when the planets are ingested. Thick regions 
on the HR diagram indicate when A?;sini> 8 km s~^. 
Of these RRPs, we find that 26 (72%) accrete planets 
before reaching the base of the RGB, and 10 (28%) do so 
while still on the main sequence. Of all the PH stars that 
ingest planets, 28 (31%) do so before arriving at the base 
of the RGB and 10 (11%) ingest planets before leaving 
the main sequence. All of the stars that ingest planets 
on the MS are transiting planets, a result consistent with 
the recent evidence that many transiting planets are in 
the midst of tidal mi gration and are not in stable orbits 
(| Jackson et al.l[2008h . Such tidally evolved planets are 
more likely to be found with transiting surveys because 
the range of inclination angles resulting in a transit in- 
creases for closer orbiting planets. 

The accretion of planets before the RGB is a signifi- 
cant result because these pre-RGB sta rs have yet to un- 
dergo rotational braking (|Gravl[l989l ). which is a mag- 
netic dynamo mechanism thought to be responsible for 
the nearly ubiquitous low rotational velocities of red gi- 
ants. In Figure I12[ this transition occurs at a temper- 
ature of ~5550 K. The proposed braking mechanism is 
triggered by the deepening convective envelope; the asso- 
ciated magnetic field forces co-rotation of escaping mass, 
allowing that mass to carry away large amounts of an- 
gular momentum. Thus, if planetary angular momen- 
tum is added to the star prior to this braking mecha- 
nism, that angular momentum will be lost before the 
star reaches the RGB. In addition, this braking mecha- 
nism should continue working throughout the RGB evo- 
lution, decreasing the amount of time that the rapid rota- 
tors spend in this high rotation state while on the RGB. 
Considering this diminishing effect of magnetic braking, 
planet ingestion probably accounts for less of the rapid 
rotator population than found in ^Ql 

Even in the absence of a magnetic dynamo, significant 
angular momentum can be carried aw a ,y b y stellar winds 
as demonstrated by ISoker fc Harpad (pOOQl , who stud- 
ied the angular momentum evolution of spun-up RGB 



stars onto the HB. The rate that stellar angular momen- 
tum is lost {J If) is proportional to the stellar angular 
velocity (o^*), the mass loss rate (Af), and R^. Ignor- 
ing interactions with companions for the moment, we 
expect the angular momentum loss to increase as the 
star ascends the RGB b ecause of mass loss. Although 
the iGirardi et al.l (|2000D evolution models we used ig- 
nore mass loss, rn odels that do include mass loss (e.g., 
ICharbonnel et aL|[T99 9) show that mass loss accelerates 
as a star ascends the RGB and is thus most substan- 
tial near the tip. Additional angular momentum is lost 
when one considers the interactions with companions. 
The tidal interactions t hat cau se planets to spiral in will 
increase , and Siess fc Livid (fl999 ) found that planet 
absorption can actually increase M on the RGB. Both 
of these will increase J. Consequently, ignoring mass 
loss will lead to increasingly inaccurate results the fur- 
ther along the stellar evolution the assumption is applied, 
which is why we chose to end our simulations at the RGB 
tip. 

To further illustrate the importance of considering 
mass loss in an analysis of planet-induced rotation in 
post-M S stars, we analyze the recent work bv lMassarottil 
(|2008f l. Contrary to our results. iMassarottil (|2008l) finds 
that HB stars have the highest probability among evolved 
stars of exhibiting detectably larger rotational velocities 
due to planet ingestion. However, this result is only pos- 
sible if the star can hold onto a substantial fraction of its 
envelope mass and any added angular momentum for the 
entire RGB evolution through the helium flash. Under 
this assumption, a star at the RGB tip can ingest plan- 
ets with relatively large Op (and thus relatively larger an- 
gular momentum) and not become a rapid rotator, but 
when that same star evolves to the HB, the stellar radius 
and moment of inertia decrease enough for that same 
angular momentum to cause rapid rotation. However, 
given that a significant fraction of the accreted angular 
momentum may be lost, the probabilities derived for the 
HB stars inay in fact be much lower than calculated in 
IMassarottil (|2008l ). Moreover, the role of planet accretion 
by RGB stars on the nature of HB stars is probably far 
more complicated than a simple calculation of rotational 
velocity from the net change in the stellar angular mo- 
mentum. The deg ree of mass loss , particularly driven by 
plane t accretion, (jPeterson et al.l[l983t ISoker fc Harpad 
120001 and references therein), and the redistribution of 
angular momentum bet ween the stellar core and enve- 
lope during helium flash (Soker fc Ha rDadl2007f ) has been 
proposed to be responsible for a rich variety of HB star 
properties, e.g., temperature, rotation, and luminosity, 
that cannot be explained by other means. 

5.3. PH Stars as a Proxy for Red Giant Populations 

In this section, we discuss the suitability of using main 
sequence PH stars as a representative sample of the pre- 
cursors to present-day red giant stars. While these PH 
stars are the obvious choice for investigating the role of 
planet accretion in spinning up giant stars, they are not 
an ideal choice. Obviously red giants are further along in 
their stellar evolution than the PH stars, which implies 
that the giants are either older on average than the PH 
stars or they are of similar age but have evolved faster 
than the PH stars. Faster evolution occurs in stars that 
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Fig. 12. — Stellar evolution tracks on the HR diagram for the RRPs. The triangles mark when the planets are ingested. Thick colored 
regions of the tracks mark rapid rotation that occurs before the RGB base in orange and from the RGB base onward in green. Panels are 
ordered by time of planet ingestion from early to late in the stellar evolution. 

ity with stellar temperature in main sequence PH stars. 
This result is evidence for primordial metal-richness of 
PH stars because hotter MS stars have shallower con- 
vective zones; if the atmospheres of PH stars were en- 
riched with metals through the planet formation pro- 
cess, then the hotter stars with shallower envelopes would 
have higher apparent metallicities than cooler stars with 



are relatively more massive or more metal poor. Among 
the global properties of the PH stars and the red gi- 
ants that may differ in a significant way, we focus on the 
metallicity, mass, and spatial distributions. 

Stellar metallicity complicates the interpretation of our 
results because of the well-defined stellar metallicity- 
plane t frequency (MPF) correlation (Fischer & Valenti' 
[200l for main sequence stars, which is the power-law 
relationship between the metallicity of a MS star and 
the likeli hood that the star has a planetary compan- 
ion. The ISozzetti et al.l ()2009 f) planet hunting survey of 
low metallicity stars extended the metallicity range of 
the correlation down to -2.0 dex. Although they con- 
firmed the power law behavior in the more metal-rich 
stars, they also found some indication that the power law 
may flatten to a constant probability at low metallicities. 
Whether this MPF correlation is "nature" or "nurture" 
is disputed. Either metal-rich stars are predisposed to 
having planets (nature) or the presence of planets can 
enrich the apparent stellar atmospheric metallicity (nur- 
ture). 

In s upport of the nature scenario, iFischer fc Valentil 
()2005[ ) found that there was no trend of stellar metallic- 



dee per envelopes. Suppor ting the nurture explanation is 
the iPasauini et all (|2007[ ) discovery that PH giant stars 
do not show a MPF correlation. The most natural ex- 
planation for the absence of the correlation is that the 
significantly deeper convective envelopes of giant stars 
have erased the planet formation pollution seen in the 
main sequence PH stars. One would expect to see the 
MPF correlation in giant stars if intrinsically metal rich 
stars are predisposed to forming planets. 

Either explanation for the MPF correlation has some 
relevance to our analysis and the interpretation of our re- 
sults. A comparison of the metallicities of the PH stars 
in FigureOwith the red giants in Figure [1] shows that the 
PH stars are on average more metal rich than the giants. 
Under the nature scenario, we consequently overestimate 
the original planet frequency in the present metal-poor 
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giant populations and overestimate the contribution of 
planet ingestion to the rapid rotation seen in these stars. 
Conversely, under the nurture scenario, the primordial 
metallicities of the PH stars are lower than those pre- 
sented in Figure [5] and more closely match those of the 
giants. However, if only the thin atmospheres of PH stars 
are polluted with metal-rich material, then we modeled 
the evolution with stellar models that were too metal- 
rich. As an experiment to see how this affects the out- 
come of our simulation, we reran our simulations, match- 
ing each PH star with stellar models that were system- 
atically more metal-poor than the measured [Fe/H]. We 
tried a range of systematic offsets between -0.1 and -2.0 
dex in increments of 0.1 dex. We find that the total num- 
ber of RRPs in our simulation varied over the range of 
36 to 42, while the average lifetime spent as a rapid rota- 
tor steadily dropped from 31.0% down to 18.4%. There- 
fore, our possible overestimation of the stellar metallicity 
because of pollution may have also resulted in an over- 
estimate of the number of expected rapid rotators from 
planet ingestion. 

The differences in the M^, distributions of the PH stars 
and red giants are also important. As illustrated in Fig- 
ure [1] the actual mass range covered by a sample of field 
red giant stars extends to higher masses than what is 
probed by the most prolific exoplanet surveys. This is 
primarily a bias in the exoplanet surveys. More mas- 
sive main sequence stars have broader spectral features 
and more atmospheric "jitter," i.e., larger amplitudes in 
the random variations of the measured radial velocities 
originating in the stellar atmosphere. Both of these ef- 
fects severely limit the attainable radial velocity precision 
that is needed to detect the period ic stellar wobbl e s that 
imply the presence of a planet. iJohnson et al.l ()2006l ) 
began a survey to circumvent these difficulties by prob- 
ing the planet frequency around more evolved massive 
stars. In particular, his survey focuses on the sub gi- 
ant branch where masses are most easily discriminated 
on a color magnitude diagram. Preliminary results from 
this survey point toward a. paucit y of close-orbiting plan- 
ets (| Johnson et al.l I2007E l200l : none of the 15 plan- 
ets discovered in this survey orbit closer than 0.8 AU. 
This observed paucity cannot be an observational bias 
because radial velocity surveys are most sensitive to the 
smallest orbital semimajor axes. Nor is this inner limit 
consistent with the maximum radius of planets that have 
tidally spiraled into the star; our predictions from Figure 
[3] suggest an inn er limit that is a t least an order of mag- 
nitude smaller. iJohnson et al.l (|2007bf ) also considered 
post-MS engulfment of planets to be unlikely. Presently, 
the most plausible explanation is that the lack of close- 
orbiting planets around more massive stars is a natu- 
ral consequence of planet formation and migration and 
gives an indication th e role of M* in these two processes 
(I Johnson et"aI]l2008^ . In our simulation, of the 50 stars 
harboring planets no closer than 0.8 AU, only 4 (8%) are 
predicted to be RRPs. Consequently, these more massive 
stars might have a value of Jrbp that is four times lower 
than that of the stars we probed. On the other hand, 
M* also plays a role in the planet frequency, and stars 
more massive than 1.3M0 are five times more likely to 
harbor planets than low-mass M dwarfs (jJohnson et al.l 
l2007al . which more than compensates for the decreased 
efficiency of spin-up. 



Finally, differences in the spatial distributions of the 
PH stars and red giants may also affect the appropriate- 
ness of comparing these two stellar samples because of 
the possibility that the locally-observed planet frequency 
is not ubiquitous throughout the Galaxy. Most known 
exoplanets were found by the radial velocity technique, 
and these surveys target nearby stars (within 100 pc), 
whereas the rapidly rotating field giants that we seek to 
explain are generally a few kpc away. One project that 
has made strides in addressing the deficiency of known 
extrasolar systems around distant stars is t he Sagittar- 
ius W indow ExoPlanet Survey (SWEEPS, ISahu et all 
l2008l ). which surveyed the Galactic bulge for one week 
using HST. The SWEEPS survey found that the fre- 
quency of planets around stars i n the bulge was con- 
sistent with those of nearby stars (|Sahu et al.ll2008l ). al- 
though the authors admit that their small sample size 
of only 16 planet candidates allows errors of factors 2-3 
in the planet frequency. Additionally, the bulge is rela- 
tively more metal rich than the solar neighborhood and 
perhaps significantly more so than many of the distant 
red giants comprising the present day rapid rotators. 

5.4. Accretion Complements Other Spin-up Mechanisms 

Angular momentum dredge-up from the stellar core is 
an oft-cited possible alternative mechanism for spinning 
up red giant stars. If giant stars are not rigid rotators 
and if their cores have high specific angular momentum, 
then the deep mixing that occurs at first dredge-up can 
spin up giants' convective envelopes. Rapid rotation in 
giant stars is only consistent with this mechanism if they 
have evolved beyond first dredge-up, whereas we have 
found that planet ingestion is most effective at creating 
rapid rotation in the earliest states of the RGB (recall 
Figures m and (12) . Unfortunately, statistics of the ob- 
served occurrence of rapid rotation as a function of RGB 
evolution time are sparse. Furthermore, any such statis- 
tics are complicated by the possibility of confusing RGB 
stars with stars in the red clump (the compressed HB 
of metal-rich stars) or AGB; all of these evolutionary 
stages overlap in an HR diagram particularly when the 
sample of giant stars have a wide range of ages, masses, 
and metallicities. Nonetheless, it is reasonable to as- 
sume that rapid rotation on the lower RGB is likely due 
to planet accretion whereas rapid rotation on the up- 
per RGB is more likely to be from angular momentum 
dredge-up because these mechanisms ar e most efficient at 
compl ementary locations on the RGB. iMassarotti et al.l 
()2008( ) found three field giants in their survey that have 
moderately high rotation, that are first ascent red giants 
that have not reached first dredge-up, and that they con- 
cluded could only be explained by planet accretion. Of 
the single giants in their sample, these three giants that 
likely accreted a planet account for 0.4% of their survey, 
and account for a slightly larger percentage of first ascent 
RGB stars. This percentage is consistent with our pre- 
diction that one would expect to find 0.58% of a sample 
of first ascent red giants to have been spun up by planet 
accretion. 

Alternatively, the disparity between the predicted 
number of rapid rotators from planet accretion and the 
larger fraction of rapid rotators seen in present-day gi- 
ants may simply be resolved by the existence of low-mass 
stellar companions. Although close-orbiting brown dwarf 
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companions are rare, close-orbit ing stellar companions 
occur in roughly 11% of stars (jGrether fc Lineweaved 
|2006[) . Consequently, some of the supposedly single red 
giant rapid rotators may in fact have undetected stel- 
lar companions with which they have tidally interacted. 
Those stellar companions on the smallest mass scales are 
the most likely to avoid detection given their relatively 
low luminosity and gravitational effects on their hosts 
compared to more massive companions. 

6. SUMMARY 

This study was motivated by the question of whether 
the phenomenon of rapidly rotating red giant stars can be 
explained solely by spin- up from the accretion of planets. 
We chose the population of main sequence planet-hosting 
stars as a proxy sample of rapid rotator progenitors. De- 
pending on the role of metallicity in the true frequency 
of planets among main sequences stars, our choice of a 
relatively metal-rich proxy biases us to deriving an up- 
per limit to the number of rapid rotators generated by 
planet absorption among the local Galactic distribution 
of red giants. On the other hand, other factors in our 
model, such as our possible overestimation of the con- 
vective envelope mass beyond first dredge-up described 
in and the choice of a relatively large value for de- 
viant rotation rates, i.e., v sin i > 10 km s^^ as deviant in 
stars for which the mean value of wsini is 2 km s~^, bias 
our results in the opposite sense and at least partially 
compensate for the other uncertainties. 

In our experiment, we used theoretical stellar evolu- 
tion tracks to follow the individual evolution of planet 
hosting stars through the first ascent of the red giant 
branch. We find that rapid rotators can be created from 
planet ingestion, though we find that the frequency of 
spin-up by planets can only account for ~10% of rapid 
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rotator giants seen in the field. This disparity may be 
resolved if some of some of the RGB stars are in fact 
AGB or red clump stars. Even if our RGB stars are 
not confused, this disparity may simply indicate that 
other mechanisms, such as angular momentum dredge- 
up or low-mass stellar companion interactions, are also at 
work. Contrary to our initial expectations, our modeling 
suggests that many of the known exoplanets will be in- 
gested by their stars prior to the RGB, particularly those 
planets in close-in orbits that have been found by tran- 
sit searches. Additionally, mass loss on the RGB as well 
as magnetic braking could significantly reduce the con- 
tribution of planet absorption to the formation of rapid 
rotators and thus require further scrutiny. Nonetheless, 
we find that the lower RGB is the best place to look 
for spin-up that is unambiguously due to planet accre- 
tion, which is possible both because the largest fraction 
of our predicted rapid rotators from planet accretion are 
rapidly rotating in this phase and because first dredge- 
up, the time of possible angular momentum dredge-up 
from the core, has not yet occurred. Finally, there is 
growing evidence that relatively massive stars are more 
likely to harbor planets than less massive stars and that 
planets around these more massive stars have a different 
distribution of semimajor axes. These trends warrant 
further study, particularly since many red giant stars are 
more massive than the planet-hosting stars used in this 
experiment. 
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